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Abstract 

In d-wave superconductors, the electronic density of states (DOS) 
induced by a vortex exhibits a divergency at low energies^ N vor t ex (E^ ~ 
l/\E\. It is the result of gap nodes in the excitations spectrum outside 
the vortex core. The heat capacity in two regimes, T 2 /T^ 3> B/B C 2 
and T 2 /T c 2 < B/B c2 , is discussed. 



1 Introduction 

Two different energy scales, Fermi energy Ep and the gap amplitude A <C 
Ep, govern the dynamics of fermions in superconductors. This leads to sev- 
eral levels of description of excitations in the inhomogeneous background of 
the order parameter produced by vortices and textures. (1) A quantum- 
mechanical approach was used for calculation of the discrete spectrum of 
bound states in the vortex core 0. This calculations revealed the existence 
of the low-energy branch 

E Q (Q) = -Quo (l.i) 



where the orbital quantum number Q is half of odd integer for conventional 
vortices. Here we consider a two-dimensional case, i.e. we neglect the depen- 
dence on the momentum p z along the vortex axis. The interlevel distance of 
bound states is small compared to the gap, typically u ~ A 2 /E F . 

(2) In the quasiclassical (or Eilenberger) approach, quantum mechanics 
is applied only for the motion along the trajectory. The energy spectrum 
is characterized by the continuous impact parameter b (or the continuous 
orbital quantum number Q = ppb) and by the trajectory direction (angle a) 
in x,y plane f2j. The low-energy part of the spectrum corresponds to the 
chiral branch in Eq.(l.l): 

E (Q,a) = -Qu Q {a) . (1.2) 

In a conventional axisymmetric vortex in s-wave superconductors, uq does 
not depend on a. In a cf-wave superconductor, the gap A(a) is anisotropic, 
it has the underlined tetragonal symmetry and exhibits four gap nodes. The 
typical example which displays these properties is A (a) oc cos(2a), which 
has four gap nodes at a$ = {2k + l)7r/4 with integer k. In the vicinity of 
each gap node A(a) ps A'\ ao (a — a ). As a result the energy spectrum in 
Eq.(1.2) has a four-fold symmetry and is given by || 

(A') 2 1 

u {a) « (a-« ) 2 i ^ln 1 , (1.3) 

tip \a — ato\ 

in the vicinity of the gap nodes. Properties of the spectrum in the whole range 
of E and a was discussed in |5[ in connection to the scanning tunneling 
microscopy experiment 0. 

The quantum limit of Eq.(l.l) is restored by using the Bohr-Sommerfeld 
quantization rule for the orbital momentum Q and the canonically conjugated 
angle a: 

/ da Q{a,E) = 27r(n + 7) , Q{a,E) = — , 1.4 

Jo ^o{ a > 

where n is an integer and 7 is of order unity. This gives the discrete levels 



E = -(n + 7) 



da 1 1 1 

2tt U1q{o) 



;i.5) 
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For axisymmetric vortices in a s-wave superconductor, Eq.(l.l) with Q = 
n + 1/2 is restored if one chooses 7 = 1/2. Actually the choice of 7 (i.e. 
either 7 = 1/2 or 7 = 0) is dictated by the symmetry of the superconducting 
state in presence of the vortex 0. 

For the d-wave case of Eq.(1.3), the integral in Eq.(1.4) diverges at angles 
close to the gap nodes. This indicates a singular behavior of the energy 
spectrum, which originates from 1/r tail of the superfluid velocity v s (r) far 
from the vortex core. 

(3) At large distances from the core, a pure classical approach can be used 
in which the coordinate r and the momentum p of the fermionic quasiparticles 
are considered as commuting variables. The energy spectrum is given by 



E = ±y t\ + A 2 (a) + k ■ v s (r) , (1.6) 

where = Vp(k — kp)- This approach has been used for calculations of 
a non-analytical behavior caused by the gap nodes. In 3 He-A, the order 
parameter texture induces the effective v s (k,r) and the point nodes lead to 
a non-analytical density of the normal component at T = ||. The same 
effect of gap nodes resulting in a nonzero DOS in presence of the conventional 
v s was discussed in P, ITU], ITTJ. Due to an inhomogeneous vortex-induced 



velocity v s (r) in the mixed state of a d-w&ve superconductor, a non- analytical 
dependence of DOS on magnetic field has been found: N(0) ~ Np(B /B^) 1 ^ 2 , 



where Np is DOS in the normal state 12, HI. The numerical factor in this 



dependence has been calculated in [13|. Possible experimental realization of 



such behavior was discussed in [14, |TJ|, Kj . 

Here we consider the energy-dependent DOS, N(E), for an isolated vortex 
and find that N(E) oc 1/\E\ at E — > using both the quasiclassical and 
classical approaches. 



2 \jE divergence of DOS in the classical ap- 
proximation. 

DOS in a homogeneous two-dimensional <i-wave superconductor is 

N(E) =2/^ \S{E T ^ + • (2.1) 
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Here the factor 2 accounts for the spin degrees of freedom and the factor 1/2 
takes care of the double counting of particles and holes with different signs 
of the square root under the 5-function. 

For small E, the main contribution comes from the four gap nodes (see 



also 1 1 



N{E) = 2 ^h lEl = 2Npl ~§ ' (2 - 2) 

where A' = d a A at the gap node and N F = m/ir is the DOS in the normal 
state. This leads to a Independence of the specific heat whose experimental 
evidences are discussed in |TJ|, |i~7fl . 

The vortex contribution comes mainly from the superflow around the 
vortex and from momenta close to the gap nodes k a ||. The locations of the 
gap nodes are not important, however we assume the tetragonal symmetry 
with k a = ±k F x, ±k F y for simplicity. All four nodes give equal contributions 
thus the extra vortex-induced DOS is 

,(E) = N F ^ J d 2 r (\E - k F x ■ v s (r)\ - \E\) . (2.3) 



N, 



d— vortex \ 



The velocity far from the vortex is v s (r) = 4>(h/2mr). Using the variable 
u = r(2mE /%kp) one has 

N d _ vortex (E) = N F ^^ J** d(J) £ du (\u-cos(f>\-u) = \n f -^- . (2.4) 

The characteristic dimension of the region near the vortex which con- 
tributes to DOS is 

r{E ) = ^- . (2.5) 

When E decreases, the size r(E) reaches the intervortex distance Rb ~ 
£(B c2 / ' B) l l 2 in the vortex lattice. For lower energies, E should be substituted 
with (hk F /2mRB), and the square-root dependence of the DOS in the lattice 
of the <i-wave vortices|| is restored. So the Eq.(2.4) holds for the energy E 
in the range A'(B/B c2 ) 1/2 < \E\ < A'. 



3 Quasiclassical approximation. 

In this approach, the radial motion is quantized. We consider only the chiral 
energy branch. The two remaining variables are: the impact parameter b (or 
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the angular momentum Q = kpb) and the angle a of the direction of linear 
momentum k in the x, y plane. DOS in this approximation is 

N d . vortex (E) = J ^9_5{E - E (Q, a)) , (3.1) 

where E (Q, a) = —Qu>o(a). Note one sign under the (^-function as compared 
to Eq. (2.1): now there is only one type of excitations on the chiral branch. 

For the d-wave case where u>o(a) is given by Eq.(1.3), the integral in 
Eq.(3.1) diverges near the nodes: (/ da/u> (a) ~ / da/ (a — a ) 2 , see 0). 
This divergence is related to the large extension of the radial wave function 
\I/(r, a) when a is close to the direction of the gap nodes. To treat this 
divergence properly, one must include \I/(r, a) into the equation for DOS 
explicitly: 

/da r 
— J d 2 r\y(r,a)\ 2 6(E-k F rsm(<f>-a)LU (a)) . (3.2) 

Here we introduced the impact parameter b = r sin(0 — a). The radial wave 
function at large distances r > ( in a vicinity of the gap node can be obtained 
from the Eilenberger equations [HI]: 

f A' \ 
\^(r, a)\ 2 = mA'\a — a | ex P — 2r — la — a | • (3.3) 

\ v F J 

Eqs. (3.2-3) solve the problem of the divergence, since they have a wider 
range of applicability than those in Ref.||, where the condition of small 
impact parameter compared to the radial extent of the wave function was 
used. Under this condition, Eqs. (3.2-3) indeed transform to Eq.(3.1). 

We integrate Eq.(3.2) over cf> introducing variables r = pvp/2A', a — a = 
a and the parameter a = LA'/\E\ (where L = — ln|a|). Summing over 4 
gap nodes, one obtains 

N d _ vortex (E) = ^ / da a pdpe~ pa —= Q(apa 2 - 1) . 

7iA'\E\ Jo Jo y(apa 2 ) 2 - 1 

(3.4) 

This integral is independent of a and is equal to tc/ 2, thus 



v 2 

N d _ vortex {E) = N f tt-^- . (3.5) 
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This is two times larger than Eq.(2.4) obtained using classical considerations. 
We can thus conclude that the classical approach is not able to give the 
correct numerical factor though it feels the relevant physics and provides a 
correct order-of-magnitude estimate. The same is known for DOS in 3 He-A 
textures: the exact DOS obtained quantum-mechanically in |l9j is two times 
the classical result of [El. 



4 Conclusion 

The singular behavior of DOS can be seen in the temperature dependence 
of the cpecific heat in presence of a vortex lattice. Being averaged over 
vortices, DOS (per one Cu0 2 superconducting layer) contains the bulk term 
Eq.(2.2) and the vortex term npNd_ VOTtex {E), where = B/§ is the flux- 
line density, and $o is the magnetic flux quantum: 

N(E) = ^-\E\ + . (4.1) 



This results in the specific heat: 

/oo 
N(E) cosh 
-oo 



00 , _ 2 / E \ E 2 dE 



IT) AT 2 




^(3)^ + 2^!^ , • (4.2) 

The second term in Eq.(4.2) is the temperature-independent linear-in-field 
correction of order NpT c B/B c2 to the dominating field-independent bulk 
term ~ NpT 2 /T c . Two terms become comparable at the temperature T/T c ~ 
{BjB c2 yi 2 , where the crossover occurs to the square-root behavior C(T) ~ 



TN F ^J B / B c2 of Ref.g at lower temperatures, T/T c < (B / B^) 1 ' 2 < 1. 
Both high-temperature and low-temperature asymptotic dependences of 



the specific heat on the magnetic field, N F T C B / B c2 and NpT^J B/B c2 , come 
due to the velocity field far from the vortex and do not depend on the details 



of the vortex core structure discussed in [[2C], [2T[ . 
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